Slow long-range decay of bound Hartree-Fock orbitals and enhancement of the 

exchange interaction and tunneling 
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Exchange interaction strongly influences the long-range behaviour of localised electron orbitals. 
It violates the oscillation theorem (creates extra nodes) and produces a power-law decay instead of 
the usual exponential decrease at large distances. For inner orbitals inside molecules decay is r~^, 
for macroscopic systems cos (fc/r)r~'^, where kf is the Fermi momentum and u — 3 for ID, 3.5 for 
2D and 4 for 3D crystal. Slow decay increases the exchange interaction between localised spins and 
the under-barrier tunneling amplitude. The under-barrier transmission coefficients in solids (e.g. 
for point contacts) become temperature-dependent. 



PACS numbers: 31.15.xr , 71.15.-m , 71.70.Gni 



I. INTRODUCTION 

The Hartree-Fock equation for an electron orbital ^'(r) 
in an atom, molecule or solid has the following form: 



2m dr^ 



■f{r) + {U{r)~E)-f{r)^K{r) (1) 



if(r) = ^*^(r) J M,^(r')t^^*(r')dr' 



(2) 



Here the summation runs over all electron orbitals ^'^(r) 
with the same spin projection as ^'(r). Now consider, 
for example, an inner electron atomic orbital Is. The 
solution of the Schroedinger equation in potential U{r) 
has a very small range Ub/Z where Z is the nuclear 
charge. Outside this range the orbital decreases expo- 
nentially as exp {—rZ/ab). In the Hartree-Fock equation 
((H) such rapid decay is impossible. Indeed, if ^'(r) ~ 
exp {—rZ/ab) the left-hand-side of eq. ([1]) would be ex- 
ponentially small while the right-hand-side is still large 
since K{y) in eq. contains higher orbitals ^'^(r) which 
have larger range. The behaviour of the inner Hartree- 
Fock orbitals inside atoms have been studied analyti- 
cally (in the semiclassical approximation) and numeri- 
cally in Ref. The dependence on the radius r can 
be found from the multipole expansion of |l/(r — r')| 
in K{t); the slowest decay normally comes from the 
dipole term r' /r"^) and/or last occupied orbital ^'^(r), 
K ^ ^'g(r) / r"^ ■ The extra nodes appear since the orbitals 
^^(r) oscillate. For example, the Is orbital in Cs atom 
has 3 nodes [l| (without the exchange term a ground 
state has no nodes) . The existence of extra nodes in so- 
lutions of Hartree-Fock equations was also mentioned in 
the book 2]. Outside the atom all orbitals decay with 
exponential factor for an external electron 

Inside solids there are electrons in the conducting band 
which occupy the whole crystal. It has been pointed out 
in Ref. [l| that this leads to a long oscillating tail of 



bound electron orbitals. The effect of the exchange in- 
teraction K{y) has been estimated in the free band elec- 
tron approximation '^q{v) = exp (iq • r). An orbital of a 
bound electron decreases at large distances as [1] 



«'(r) 



;(/c/r)/r 



(3) 



where kf in the Fermi momentum. 

The derivation of this expression assumes the presence 
of a partly filled conducting electron band. However, in 
atoms and molecules of any length the exchange enhance- 
ment of the inner orbital tail may be mediated by a com- 
plete electron shell. The question is: can the exchange 
enhancement in solids be mediated by a nonconducting 
electron band? A special interest in this problem may be 
motivated by spintronics and solid state quantum com- 
puters based on spin qubits. The long-range tail of the 
wave function could, in principle, lead to an enhance- 
ment of the exchange spin-spin interaction between the 
distant localised spins, and enhancement of the under- 
barrier tunneling amplitude. 

A special feature of the "long-tail" mechanism is that 
the state of the band electrons does not change, i.e. there 
is no need to have polarization of the conducting band by 
the localised spin. The mediating band electrons produce 
the mean exchange field K{t) in eq. ^ only. Therefore, 
this "long-tail" effect is different from other effects like 
the RKKI interaction [4j| and the double exchange spin- 
spin interaction suggested by Zener Q (see also Refs. 
[a, 0| and description of Anderson and Kondo problems 
, e.g., in the book 

To investigate this problem in the present paper we 
perform calculation of the tail using the Bloch waves and 
tight-binding band electron wave functions. 



II. ATOM 

Let us first explain how the long tail appears in atoms 
1]. The radial equation for a Hartree-Fock electronic 
orbital ^i(r) = r(j)i[r) is 



2m dr^ 



(4) 



2 



U 



(5) 



The radial exchange term can be obtained using the mul- 
tipole expansion of l/|r — r'|. Outside the radius of an 
inner orbital S,i (e.g. in the area r > ub/Z for Is) 



E 

k>0,n 



(6) 



dinner (r) « coust ^valence (r)/r2 . In macroscopic systems 
there is a large number of electrons occupying the valence 
band and the contribution of different valence electrons 
interfere in the exchange term in Eq.Q. This interfer- 
ence changes the long range behaviour. 

The equation for a bound electron wave function ^'^(r) 
in a crystal contains the exchange term from Eq.([2) de- 
scribing the exchange interaction of the bound electron 
with 2F mobile electrons: 



Here Cnk are the standard angular momentum dependent ^ ^2 ^2 
coefHcients and 6„fe /r'''^„(r)^i(r)(ir. For the multipo- K{r) ^ g{r - r')[- ]5'(,(r')dr', (10) 

laritv k = thp intpcrra] /)„ i. = f) diip tn tVip nrthncnnalil^v |r r | 7" 



larity fc = the integral bnk = due to the orthogonality 
of radial wave functions with the same angular momen- 
tum. 

Now we can discuss the large distance behaviour of the 
orbital £,i{r). We will use Is orbital in Xe atom {Z = 54) 
as an example. The last occupied shells are ...5s^5p^. 
The orbital 5s does not contribute to Ki{r) since in this 
case the multipolarity of the exchange integral is fc = 
and the orthogonality condition makes bnk = 0. The 
exchange integral ls5p has fc = 1, therefore, at r ^ as 
and outside the atom Kis{r) « C'5p.i65p_i^5p(r)/r^. 

The solution of Eq. ^ may be presented as 



e,(r)=4/'--(r)+er(r) 



(7) 



Outside the radius of the inner orbital (r > ag /Z for 
Is) the energy Ei is much larger than other terms in the 
denominator of Eq. ([5]) which are of the order of En 
(since the opertor in the denominator acts on ^„). In our 
example the energy of Is is \Ei\ = Z'^ x 13.6 eV=4 • 10'' 
eV while the 5p energy is \En\ ~ 10 eV. Therefore, we 
can approximately write 



er(r) 



K^{r) 



U, 



eff 



E, 2m{Uef f~Ei) dr^ Uef / - E, 



(9) 

The free solution in this area may be described 
by the semiclassical (WKB) approximation, 
^^^'^^(^r) ^ IpI^^/"^ exp{~ J \p\dr/h); it has the usual 
range as/Z = 0.02aB for Is. Comparison with the 
numerical solution of the Hartree-Fock equation for Is 
orbital has shown that within ~1% accuracy it is enough 
to keep the first two terms in the expansion Eq. ([9|) 
beyond the classical turning point, and only one term 
at r > lOaB/Z. Similar results have been obtained 
for the Dirac-Hartree-Fock orbitals which include the 
spin-orbit interaction and other single-particle relativis- 
tic corrections [l|. Thus we see that at large distances 
£_is{r) « const ^5p(r)/r2. 



III. ID, 2D AND 3D SYSTEMS 

If we consider a molecule instead of atom, in- 
ner electron orbital will behave the same way. 



ff(r-r')^^vl/„(r)M/„(r')t. 



(11) 



Summation goes over F mobile electron states ^'n(r) with 
the same spin projection. To account for the orthogo- 
nality condition / \I'„(r')t\I/f,(r')(ir' = in Eq. we 

excluded the zero multipolarity term from the Coulomb 

2 22 
integrals, replacing ^r^rpj by ^rzrpj — In the "exact" 

expression (fTU]) the subtracted term disappears after 

the integration over r' since J ^'„(r')^^'I'b(r')(ir' — 0. 

Let us start discussion of crystals from the simplest 
problem - a ID chain of N atoms separated by distance 
a. The wave function of a mobile electron can presented 
as 



4'„(r) = L-i/2e^'="^«fc(r), 



(12) 



where Vk (r) is a periodic function in x-direction and L = 
Na is the length of the chain. To perform the summation 
in Eq. (|lip analytically we neglect dependence on fc in 
Vk (r) . Taking the standard set of the wave vectors fc„ — 
27rn/L, n = 0, ±1, ±q, where F — 2q+ 1, we obtain 

/ /X / X / ,x sin [fcf (a; — a;')l , „x 

g{r - v') = v{v)v{r') , . (13) 

X — X 

where kp = fir /a and f — F/N is the band filling fac- 
tor. Now we can find the exchange term Eq (fTO|) . The 
leading term in the multipole expansion (r' << r) of 



|r— r'l r r- 

for K{r) at large distance: 



leads to the dipole approximation 



K{r) 



sMk,x)Jx'cosik,x'Hr')MrVr' 
los {kfx) / x'sin(fc/a;')w(r')*6(r')'^r'jl4) 



It is easy to extend the problem to 2D and 3-D cases. 
In 2D case we obtain 



, , Ji{kfR) sin(fc/i?-7r/4) 
g{v-v)^v{v)v{v)^J^ ^^3^, (15) 

where R = r — r' and Ji is the Bessel function. In 3D 
case 

t;(r)i;(r')r /. x sin(fcfi?), , , 



27r2i?2 



kfR 



3 



Substituting these results into Eq ([TO]) we obtain in the 
dipole approximation that the exchange interaction term 
decays as 



K{r) ^ cos {kfr)r 



where kf is the Fermi momentum and ly = 3 for ID, 
1/ =3.5 for 2D and u =4 for 3D crystal, i.e. v = {5 + d)/2 
where d = 1,2,3 is the dimension. 

Note that the expressions (|14ll7p do not vanish if the 
electron band is complete. Instead they have fast oscil- 
lations if the electron Fermi momentum fc/ is large. This 
conclusion looks surprising since a complete band does 
not contribute to the conductivity. If this conclusion is 
correct, one may have an enhanced tunneling amplitude 
or enhanced exchange interaction between distant spins 
(power suppression r"'' instead of exponential suppres- 
sion) even in non-conducting materials. However, this 
phenomenon may be an artefact of the Hartree-Fock ap- 
proximation for the wave function (antisymmetric prod- 
uct of the Bloch waves). We, in fact, have assumed that 
the electrons in complete band are not localised, i.e. they 
are still described by the Bloch waves (fT2|) spread over 
all the crystal. Electron correlations may kill this effect, 
for example, by creating the Mott insulator where all 
electrons are localised. On the other hand, it would be 
incorrect to say that the long-range exchange is impossi- 
ble in principle. For example, correlation corrections do 
not prevent valence electrons from being present on all 
atoms in a molecule where there is no conductivity. 

The long-tail effect does not appear in any approach 
where the exchange interaction is replaced by an effective 
potential or by a density-dependent potential, e.g. in the 
density functional approach. Approximate calculations 
may also lead to other incorrect conclusions. For exam- 
ple, the long-range tail for a complete band case does not 
appear in the tight-binding approximation for the elec- 
tron wave functions. In the tight-binding approximation 
a wave function of mobile electron is 



(18) 



stitution of V&n from Eq. p8|) into Eq. (jlip and summa- 
tion over n gives the following results: 



(17) 5(r-r') = ^S(F,;-m)«'i(r-Za)4'i(r'-ma)^ (19) 



BiFJ) 



exp {i2TrlF/N) - 1 
N{exp{i27rl/N) - 1) 



exp(i7r/Z) . 
sm tt/Z , 

TTl 



where 'i'i(r — la) is the one-site wave function. The sub- 



(20) 

where I > 0, f = F/N is the band filling factor and the 
last expression is obtained for / N. For I = we have 
B{F, 0) /. Substitution of ^(r - r') from Eq. 
into Eq. (fTU]) shows that if the band is partly filled, the 
tight-binding approximation leads to the same conclusion 
K{r) ~ cos (fc/r)r^''. However, for the completely filled 
band / = 1 and sin {nfl) = 0. This means that the long- 
range exchange term vanishes in the absence of mobile 
carriers, electrons or holes. The explanation is simple: in 
the tight-binding approximation the complete band wave 
function made of the running waves Eq. (jl8p is equal 
to the antisymmetrised product of the localised electron 
wave functions 'I'i(r — ^a). The exchange interaction with 
the localised electrons does not produce the long-range 
tail. To compare with the Bloch wave expression one may 
say that the tight-binding result for the complete band 
corresponds to K{r) ^ sin (kfr) = for r = la. However, 
the oscillations of K{r) do not lead to vanishing of its 
effect on the wave functions - compare with the solution 
for atomic orbitals in the previous section. 

At finite temperature conducting electrons and holes 
appear. This activates the long-tail mechanism even in 
the tight-binding approximation and makes the under- 
barrier transmission coefficient temperature dependent. 
Here it may be appropriate to recall that a temperature 
dependence of the transmission coefficient has been ob- 
served near the "0.7 (2e^//i) structure" in the point con- 
tact conductance measurements [§, [l3| . 
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